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ABSTRACT 

The asymptotic giant branch (AGB) of the globular cluster NGC 4372 appears to 
extend to unexpectedly high luminosities. We show, on the basis of proper motions 
and spatial distribution, that the extended AGB is indeed a likely part of the cluster. 
We also present the first spectra of the very cool (2600 K), very luminous (8000 L ), 
very dusty, oxygen-rich, purported long-period variable stars VI and V2 that define 
the AGB tip. In particular, on the basis of their radial velocities, we conclude that 
VI and V2 are probably members. We find that VI and V2 are likely undergoing the 
superwind phase that terminates their nuclear-burning evolution. We hypothesise that 
the mass-loss processes that terminate the AGB are inhibited in NGC 4372 due to a 
lack of atmospheric pulsation and the high gas-to-dust ratio in the ejecta, leading to a 
delay in the associated enhanced mass loss and dust production. Previously predicted, 
but never observed, this explains the high mass of the white dwarf in Pease 1 in 
Ml 5 without the need to invoke a stellar merger. If commonplace, this phenomenon 
has implications for the mass return from stars, the production of carbon stars and 
supernovae through the Universe's history, and the AGB contribution to light from 
unresolved metal-poor populations. 
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1 INTRODUCTION 

The asymptotic giant branch (AGB) of globular clusters 
rarely extends more than a magnitude beyond the red 
giant branch (RGB) tip, either in the near-infrared or 
in bo lometric luminosity (|Bover et alj 120091 : IValenti et al.l 
2010). This is particularly tr ue in low-metallicity clus- 
ters (e.g iMcDonald et ail l201Ch . which are typically older 



|Marfn-Franch et al. 20091 ) . with less-massive stars which 



have lower AGB-tip luminosities (|Marigo et al.l l2008). 

NGC 4372 is one of the most metal-poor Galactic glob- 
ular clusters ([Fe/H] = -2.17), and lies at 5.8 kpc distance 
behind moderate optical obscuration (E(B — V) — 0.39 mag) 
at I = 301°, b = -10° (|Harrisl Il996h. At an estimated age 
close to that of 47 Tuc and M15 l|de Angeli et aLlboOSl ), its 
asymptotic giant branch (AGB) star s initially had a mass of 
«0.89 M (|McDonald et al.ll2011bh . placing the AGB and 
RGB tips at similar luminosity. Nevertheless, there are sev- 
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eral stars toward the cluster which appear to form an AGB 
extending ~3 K s magnitudes above the RGB tip (Figure 
[p. VI and V2 are variable stars in the cluster's vicinity 
( Clement e t al. 2001), and appear to form the AGB tip, but 
nothing appears published about these stars except photom- 
etry from all-sky surveys. We herein seek to confirm both 
their membership and the presence of an extended AGB. 



2 VI AND V2 
2.1 Observations 

To determine the nature of VI and V2, we took optical 
spectra using the grating spectrograph on the 1.9-m Rad- 
cliffe Telescope at SAAO. The spectrum of VI was taken 
on 2012 Apr 17 and covers 643-928 nm at a resolution of 
R — A A/ A w 4800. The low-resolution spectrum of V2 was 
taken on 2012 Apr 16, covering 608-998 nm at R « 3600; 
the high-resolution spectrum on 2012 Apr 13, covering 832- 
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Figure 1. 2MASS (J-K S ),K S diagram for stars <12.5' from the 
centre of NGC 4372 (red points). For comparison, u) Cen (cyan 
points; [Fe/H] -1.6) and 47 Tuc (green points; [ Fe/H] = -0.7 ) 
are also shown, reddened and scaled appropriately <HarriJl996l) . 
The parallelograms denote regions chosen to represent the clus- 
ter's RGB and AGB. The black line show s a [Z/H] = -2.17, 11 
Gyr Padova isochrone (Marigo ct al. 2008), reddened and shifted 
to match the cluster. The large magenta circle shows the AGB 
outlier mentioned in Section [3] 

894 nm at R « 24 000. We show these spectra along with 
literature photometry in the top panel of Figure [2] 

Both stars are clearly oxygen rich and very cool, show- 
ing TiO and VO bands: we estimate the spectral type of VI 
and V2 to be M7-M8 and M8-M 9, based on the Fluks spec- 
tral standards l|Fluks et al.lll994T ). These bands obscure the 
atomic lines, making accurate radial velocity determination 
difficult. The heliocentric radial velocity for V2 is found to 
be v T = 100 ± 10 km s _1 based on comparison to atmo- 
spheric lines, the CuNe arc lamp spectrum and subsequent 
observation of the M0 star HIP 56619. The low-resolution 
of the VI spectrum precluded an accurate radial velocity 
measurement: we estimate it lies between v r = 50 and 100 
km s _1 . These velocities are significantly different from that 
of the contaminating Galactic foreground population, which 
we measure to be 3 ± 21 km s _1 (standard deviation) on 
the basis of an initial reduction of archival VLT/FLAMES 
data of the cluster's RGB stars. 

As no variability information (and precious little reli- 
able optical data) is available in the literature, we obtained 
archival imagery taken with the ESO 2.2-m/WFI at La Silla. 
This contains t/BWc-band data with four epochs between 
2001.157-2002.466. The variability between epochs was ob- 
tained by comparison to nearby stars. We find the variabil- 
ity of V2 to be AB, V, I c > 0.85, 0.81 and 0.46 mag, re- 
spectively. For VI, we find that AB, V,I C > 0.67, 0.49 and 
0.23 mag, respectively. There is insufficient data to constrain 
any period, which may be longer than the temporal cover- 
age. Most likely, these values have been underestimated by a 
factor 1.5, but our lower limits are (respectively) thrice and 
twice th at found in the slightly -lower metallicity star K825 
in M15 jMcDonald et al.ll2010T l. 

We show our estimated optical magnitudes along with 
literature infrared data in the lower panel of Figure [2] It 
is clear that the infrared photometry does not follow a 
Rayleigh-Jeans tail, but exhibits an infrared excess. We 
have attempted to model this excess usi ng the dusty ra- 
diative transfer code |Nenkova et al.lll999h for metallic iron, 
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Figure 2. Top panel: spectra of VI and V2 from the SAAO. 
Bottom panel: S EP of V2. The hollo w points show literature 
data from DE NIS iEpchtein et al.ll994l) . 2 MASS iSkrutskie et al] 
l200rJ) . WI SE llCutri et alj|2012h . AKARI jlshihara et al.ll2010l) . 
and IRAS l lMoshir et al .119901) . The solid points show magnitudes 
estimated from archival ESO 2.2-m/WFI observations. The black 
dotted line is a 2600 K, log(g) = 0, [Fe/H] = -2 BT-Settl model 
atmosphere, reddened using E(B — V) = 0.39 mag. The inset 
shows DUSTY model fits for metallic iron (magenta, dotted line), 
amorphous carbon (blue, dashed line), large silicate grains (red, 
solid line). 

amorphous carbon and large silicate grains by assuming a 
radiatively-driven wind (Figure [2] bottom panel). It is clear 
that any of the three dust species can be made to model a 
featureless modified blackbody, though we cannot rule out 
(without 10- /xm spectroscopy) the presence of a small silicate 
feature in either emission or absorption. Assuming a lumi- 
nosity of 8000 Lq, all three models require mass- loss rates 
of M>10' 5 ^(ip/30 000) M yr~\ or expansion velocities 

of ?Joo<0.5-y/ (30 000/f/>) km s~ , where ifi is the gas-to-dust 
ratio, which we take to be 30 000 at the metallicity of the 
cluster. While these values m ay not be accurate in a n ab- 
solute sense (see discussion in lMcDonald et aPl2011ah . the 
mass- loss rate appears to be ~10x larger than the highest 
mass-loss rates in other clusters, which have be en computed 
using similar methods (jMcDonald et al.ll2011al ldf). 

2.2 Cluster members? 

The late spectral types of VI and V2 force them to be 
foreground dwarfs, optically-enshrouded young stellar ob- 
jects (YSOs), or giant stars. Their radial velocity, variabil- 
ity, and lack of strong 7665 and 7699 A K I lines preclude 
either dwarfs or YSOs. While the dark nebula Sandqvist 
149 (the 'Dark Doodad'; ISandqvistlll977< ) lies near the clus- 
ter, it is not actively star-forming nor gives extra extinc- 
tion to these stars compared to the rest of the cluster. Its 
distance is unknown, but is likely to be <500 pc, based 
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Figure 3. Reduced \ f° r fits between the spectrum of V2 and 
BT-Settl model spectra at [Fe/H] = -2. The colour scale shows 
the x 2 goodness-of-fit with darker colours being better matches. 



on a Galactic molecular gas scale height of 81 pc 
2005) . The stellar pulsatio ns suggest a luminosity of >680 
L0 ( McDonald et alll2012h . placing their distance from the 
Sun at >1.5 kpc, >260 pc below the Galactic Plane. This 
supports the idea that they are distant giants lying out of 
the Galactic Plane. 

VI and V2 are extra-ordin arily cool for stars which we 
suppose are at [Fe/H] = -2.17 (|Harrislll99fj l. Strong molec- 
ular bands preclude us from using the normal indicators 
to determine their metallicity: the Call, Fel and Til lines 
normally visible in cool star spectra are practically indistin- 
guishable from the molecular absorption, though this may 
be expected from a metal-poor star. This makes it difficult 
to use conventional techniques to estimate the stellar pa- 
rameters, thus we use three different ad hoc methods. 

Firstly, we simply note that R Leo provides a close 
match to both the detailed spectra and broad SED of V2. 
R Leo has a luminosity of 8200 L and is at a distance o f 
110 ± 9 pc (|Fedele et al.ll2005l ; Ir omano fc MatteuccHl2007t) . 
Scaling R Leo to match the near-IR flux of V2 implies a 
distance of 5 ± w 1 kpc (cf. NGC 4372 at 5.8 kpc). 

Secondly, we compare the stellar SEDs to BT- 
Settl and marc s model spectra (Gustafsson et alj 1 19751 ; 
lAllard et all 120031: iGustafsson et alT T2008) using the code 
of iMcDonald et al l ^Ooi T Assuming both stars are at 5.8 
kpc, we find that Lyi, V2 ~ 7000, 8000 Lq. These values 
are subject to a substantial error as the derived stellar tem- 
peratures (Tvi, v2 ~ 2300, 2000 K) are outside the range of 
validity of the models (2600 and 2500 K, respectively). Sim- 
ple trapezoidal integration of the SED, however, suggests 
the luminosities are no more than 15 per cent in error. 

Finally, we perform a \ 2 fit between the normalised 
observed and model spectra in the range 8470-8930 A. A 
strong degeneracy exists between temperature and surface 
gravity, and we are again limited by the range of the models, 
thus we cannot fully constrain the parameters. Furthermore, 
neither model set correctly reproduces the observed VO 
bands. For V2, the best-fitting BT-Settl model at [Fe/H] 
= -2 (Figure El is at the lowest temperature (T = 2600 K) 
with \og(g) — 0.5. The (incomplete) grid of MARCS mod- 
els shows a poorer fit at 2500 K than 2600 K. Based on 
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Figure 4. Source density within a given radius of NGC 4372, 
showing AGB stars (red line), RGB stars (grey line, scaled down 
by a factor of 30) and foreground control field (blue line, scaled 
down by a factor of 10), using the boundaries defined in Figure 
[T] Vertical lines show Poissonian errors on the number of stars. 
VI and V2 are indicated by arrows, though they lie at opposite 
sides of the cluster. 



these values, we therefore estimate the parameters of V2 as 
L = 8000 L 0J T = 2550 K, R = 460 R Q , log(p) = -1.1 
and v esc = 23 km s _1 . The best-fit match for V2 gives 
L = 7000 L , T = 2700 K, R = 380 R , log(p) = -0.9 and 
v C sc = 26 km s _1 . All these values are subject to substantial 
but presently-unquantifiable errors, but are consistent with 
parameters expected for AGB stars of that luminosity. 

However, if VI and V2 are members, we are left with 
a slight puzzle: how to reconcile the radial velocity of V2 
(100 ± 10 km s" 1 ) with that of the cluster (7 2.3 ± 1.2 km 



_|Harrislll996l; dispersion 4.6 ± 0.5 km s MGeisler et al.l 
1 1995 ). One possible solution is radial pulsation: for exam- 
ple, the change in observed radial velocity caused by pul- 
sation e xceeds ±10 km s" 1 for the strongest variables in 
47 Tuc dLebzelter et al.ll2005l) and is ±13.5 km s _1 for the 
CO bands of R Leo l|Hinklelll978l ). This could bring the 
discrepency within the error budget. Binarity might also ex- 
plain the difference, th ough globular cluste r binary fractions 
are typically very low (jlvanova et al.ll2005l ). 



3 THE EXTENDED AGB 

Using the definitions in Figure [TJ the extended AGB within 
12.5' is comprised of seven stars (listed in Table [TJ plus one 
outlier which we do not consider further. The optical and 
infrared colours are consistent with a sequence of progres- 
sively redder giant stars. This sequence starts, as expected, 
from the RGB tip. The Padova isochrone in Figure [T] shows 
the RGB tip to lie at K s « 8.4 mag. The isochrone also fol- 
lows stars through the thermal-pulsating AGB (TP- AGB), 
but does not correctly model the proposed upper AGB. 
We co uld not fit any other isochrone from iMarigo et al.l 
(2008) to this sequence. Stars above sa Q.9 Mp should become 
carbon-rich at this metallicity (see also lKarakas fc Lattanziol 
l2007l:lMarigo fc Girardll2007h . Since VI and V2 are oxygen- 
rich, we cannot expect younger (higher-mass) metal-poor 
isochrones to correctly reproduce a sequence of oxygen-rich 
stars. 

The proposed upper AGB, as defined in Figure [T] shares 
a radial density distribution with the RGB, but not the fore- 
ground population (Figure [4| . In combination with the ra- 
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Table 1. Data for the candi date upper AGB st ars. Pho- 
tometry comes from APASS llHenden et ail I2012T I; DENIS 
Epchtein et aljl994h ; 2MASS jSkrutskie et alj2006h : and WISE 



Cutri et ail l2012t) . Proper motions are shown from PPMXL 
Roeser et ail l2010h; UCA C4 (Zacharias et ai, in prep.) and 



SPM4 jGirard et alJlioil 



2MASS Photometry (mag) Proper motion (mas/yr) 

RA & Dec B V R I PPMXL 

(J2000.0) J H K UCAC4 

Radius (') [3.35] [4.60] [11.6] [22.1] SPM4 



12 28 05.686 
-72 45 19.94 
11.92 (V2) 



16.040 
7.136 
5.610 



15.470 
6.125 
4.980 



13.465 
5.483 
4.010 



11.455 
3.180 



0.7± 
2.5 ± 
-5.7 ± 



11.8 
4.0 

2.8 



12 23 00.611 
-72 40 03.04 
12.29 (VI) 



17.085 
7.383 
5.545 



15.248 
6.307 
5.094 



14.490 
5.794 
4.520 



10.227 
3.945 



3.9 ± 
18.3 ± 
14.6 ± 



12.3 

4.1 

1.6 



12 25 22.538 
-72 33 17.48 
6.48 (#3) 



15.030 
7.771 
6.207 



13.052 
6.744 
6.349 



12.195 
6.340 
6.140 



9.871 
5.984 



-17.8 ± 
10.9 ± 
■ ■ -± 



10.1 
2.9 



12 24 57.740 
-72 39 04.98 
3.58 (#4) 



14.144 
7.866 
6.457 



12.096 
6.917 
6.572 



11.273 
6.557 
6.418 



9.831 
6.258 



5.5 ± 
5.3 ± 
0.4 ± 



9.1 
1.6 
1.3 



12 25 27.364 
-72 38 41.41 
1.60 (#5) 



13.740 
8.333 
6.987 



11.877 
7.416 
7.069 



11.132 
7.131 
6.979 



9.913 
6.826 



-6.6 ± 
-9.5 ± 
-5.6 ± 



11.8 

1.6 

1.3 



12 27 19.786 
-72 35 17.53 
8.23 (#6) 



14.157 
8.858 
7.614 



12.374 
8.017 
7.703 



11.662 
7.750 
7.587 



10.399 
7.445 



-4.7 ± 
3.0 ± 
-1.9 ± 



10.1 

1.4 

1.0 



12 24 49.536 
-72 43 00.74 
5.41 (#7) 



14.081 
9.036 
7.778 



12.275 
8.203 
7.771 



11.580 
7.956 
7.644 



10.355 
7.454 



-2.6 ± 
-5.3 ± 
-5.0 ± 



11.6 

1.5 

1.1 



-3.1 ±11.8 
-4.5 ±4.0 
2.1±3.1 



-19.9 ±12.3 
-9.1 ±4.0 
12.4±1.7 



5.2 ±10.1 
-0.8 ±1.8 
...±... 



24.6 ±9.1 
12.1 ±1.6 
16.0±1.3 



-26.5 ±11.8 
-5.1 ±10.6 
-4.6 ±1.4 



3.0 ±10.1 
-3.8 ±1.5 
-0.9 ±1.0 



-20.5 ±11.6 
-7.9 ±1.9 
4.9 ±1.2 



dial velocity limits placed on VI and V2, this argues that 
these stars are indeed part of the cluster. 

Proper motions of cluster members should typically be 
<1 m as yr _1 from the cluster mean (cf. Ivan Leeuwen et ail 
2000). While the proper motions listed in Table [T] appear 
greater, the formal errors are likely under-estimated in this 
crowded field: the RGB population in UCAC4 and SPM4, 
for example, has a mean and standard deviation of Aa = 
-2 ± 30, AS = 3 ± 31 and Aa = -3 ± 21, AS = 2 ± 20 mas 
yr~\ respectively. Only star #4 shows a non-zero proper 
motion that is consistently among the catalogues, and even 
that is well within the standard deviations listed above. 



4 DISCUSSION 



Of th e 13 Galactic globular clusters with [Fe/H] < -2 (|Harrisl 
1996), only NGC 4372 has a visually obvious extended 
AGB. Many are sparsely populated enough that an extended 
AGB would be a stochastic phenomenon, with only one or 
two stars. Other extended AGBs may therefore have been 
missed. However, we note that the most-metal- poor glob- 
ular cluster (M15; [Fe/H] = -2.37; lHarrislll996l ) unusually 
contains a planetar y nebula, Pease 1, with a white dwarf 
of 0.58-0.62 M jBianchi et ail 1 19951 ; lAlves et ail |2000| ; 
iBianchi et aill200 ll ). Thi s is rather higher th an the 0.53 M Q 
found in other clusters (|Kalirai et al 1 12009). but co nsistent 
with the core mass estimated for NGC 4372 V2 (|Blockerl 
1993). Some stars may therefore reach unexpectedly-high 



luminosities in M15 too. However, NGC 4372 may simply 
be unique. 

We consider four ways an extended AGB could oc- 
cur: (1) a young, presumabl y metal-rich su b-population 
with higher mass (cf. w Cen; iLee et ail Il999l ); (2) a sub- 
population arising from evolved stellar mergers; (3) excess 
luminosity above the canonical core -mass-luminos ity rela- 
tion (cf. hot bottom burning; HBB; lHerwi3l2005l ); (4) re- 
duced mass loss during the star's evolution. 

There are four stars on the AGB between 8.0 < K s < 
7.0 mag and 26 stars on the RGB between 9.4 < K B < 8.4 
mag, giving a drop in source density of a factor of 6.5 ± 3.5 
across the RGB tip. Other magnitude ranges give similar 
results. While Poisson noise is significant, this is commen- 
surate with the factor of ~4 drop expected if all RGB stars 
survive onto the upper AGB (cf. ISandauist fc Boltei liool ) 
and suggests that population giving rise to the extended 
AGB comprises >40 per cent of the cluster. While the main 
sequence is broadened, t his is attributable to differential red- 
dening over the cluster (jAlcaino et al.lll99ll ). rather than an 
age or metallicity spread within it. To yield 0.62 Mq white 
dwarfs, such a population would also have to be exception- 
ally massive (~2 m©) and very young (~1 Gyr) to reach 
these luminosities (Blo ckerll 19931 ; ?). This negates case (1). 

The drop in sources over the RGB tip also contradicts 
case (2) , which can be rejected due to the cluster's low den- 
sity. NGC 4372 has one of the lowest central concentrations 
of stars of any globular cluster and i s not core co llapsed, 
resulting in fewer stellar encounters (|Harrisl fl996l ). While 
NGC 4372 does show an extended blue horizontal branch 
(HB) or blue straggler population (signs of stellar mergers), 
these st ars are very few in number compared to the total 
cluster (|Kaluznv fc Krzem inski 1993). In order for >40 per 
cent of the cluster to reach this phase, at least four of every 
seven stars must merge. We therefore rule out case (2). 

The temperatures requir ed for HBB should not o ccur 
at all in stars below <1.5 M© (| Ventura fc Marigd[20ir3 ) . We 
therefore do not consider case (3) a likely option. 

We are therefore left with case (4). The conventional 
mechanisms by which oxygen-rich A GB stars lose mass, al- 
ready strained at solar metallicity (Woitkc 2006]) become 
even more difficult at low metallicity: the pulsations levitat- 
ing material from the stellar surface are typically weaker, the 
dust that condenses in the outflow forms a smaller fraction of 
the wind, meaning that stellar radiation pressure has more 
difficulty accelerating both it and (by collisional coupling) 
the surrounding gas from the star. It has been suggested 
l|Bowen fc Willsorj H99ll ; IZiilstral liool ) that, at low metal- 
licity, the inability to effectively drive a dust-driven wind is 
sufficient to decrease mass-loss rates substantially. Stars are 
then given more time to ascend the asymptotic giant branch 
(AGB) before an enhanced superwind ejects the remainder 

of the envelope into the interstellar medium (ISM). 

The blue nature of the HB l|Kaluznv fc Krzeminskil 

Il993l ) suggests that pre-AGB mas s loss occurs normally . 
The remaining envelope (wO.l Mm ; iMcDonald et aill2011cl ) 
should mostly be lost via pulsation-enha nced, dust-driven 
winds at lu minosities above ~700 Lq (|McDonald et all 
l2011dl . [20121 ). However, we have seen that dust is not pro- 
duced in great quantity by any stars except VI and V2 (Ta- 
ble [IJ, and VI and V2 appear to be the only long-period 
variable pulsators (N. Matsunaga, private communication). 
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We therefore theorise that the lack of pulsation and dust 
formation delays the usual dust-driven wind until the stars 
reach much cooler temperatures. 

If this is typical of very-low-metallicity populations, this 
would have significant consequences. The increased AGB 
lifetime would lead to an increased white dwarf mass and 
decreased return of metal-rich material to the ISM. That 
similar extended AGBs are not seen in other clusters of sim- 
ilar metallicity may point to detailed elemental abundances 
being more important in iron content. 



5 CONCLUSIONS 

We provide evidence that mass loss in NGC 4372 is delayed 
due its low metallicity, explaining the extension of its AGB 
to «8000 Lq. We present spectra of the two unusually cool 
and dusty stars, VI and V2, which form the AGB tip. We 
also give evidence to associate the extended AGB, and VI 
and V2 in particular, with the cluster. While NGC 4372 thus 
far remains unique, we reconcile its extended AGB with the 
existence of the planetary nebula within M15. We finish by 
discussing the considerable impact an extended AGB has 
on stellar evolution, our view of the early Universe and its 
chemical evolution to date. 
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